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ABSTRACT 

We present a new method of measuring the power spectrum of initial perturbations to an 
unprecedently small scale of ^ lOh^^ kpc. We apply this method to a sample of 4500 Ly-a 
absorbers and recover the cold dark matter (CDM) like power spectrum at scales > 300/i^^kpc 
with a precision of ~ 10%. However at scales ~ 10 — 300/i^^kpc the measured and CDM-likc 
spectra are noticeably different. This result suggests a complex inflation with generation of 
excess power at small scales. The magnitude and reliability of these deviations depend upon 
the possible incompleteness of our sample and poorly understood process of formation of weak 
absorbers. Confirmation of the CDM-like shape of the initial power spectrum or detection of its 
distortions at small scales are equally important for widely discussed problems of physics of the 
early Universe, galaxy formation, and reheating of the Universe. 

We use the Zel'dovich theory of gravitational instability to derive statistical description of the 
properties of observed structure. Our method links the observed mass function of absorbers with 
the correlation function of the initial velocity field and therefore it avoids the Nyquist restrictions 
limiting the investigations based on the smoothed fiux or density fields. This approach is in 
general consistent with numerical simulations of the process of structure formation, describes 
reasonably well the Large Scale Structure observed in the galaxy distribution at small redshifts 
and emphasizes the generic similarity of galaxies and absorbers. 

The physical model of absorbers adopted here asserts that they are formed in the course of 
both linear and nonlinear adiabatic or shock compression of dark matter (DM) and gaseous mat- 
ter. It allows us to link the column density and overdensity of DM and gaseous components with 
observed characteristics of absorbers such as the column density of neutral hydrogen, redshifts 
and Doppler parameter. At scales > l/i~^Mpc all characteristics of the DM component and, in 
particular, their redshift distribution are found to be consistent with theoretical expectations for 
Gaussian initial perturbations with a CDM-like power spectrum. 

Subject headings: Cosmology: observations - large scale structure - qusars: absorption lines 
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1. Introduction 

One of the most important problems of mod- 
ern cosmology is to determine the power spec- 
trum of primordial perturbations. This problem 
is closely connected with physics of the early Uni- 
verse, properties of the inflation process, reion- 
ization of the Universe at moderate redshifts and 
formation of the Large Scale Structure (LSS) ob- 
served in deep galaxy surveys among others. 

Presently various observational approaches are 
used to measure the power spectrum of the ini- 
tial density perturbations. The amplitude and 
the shape of initial power spectrum on scales 
> 10/i~^Mpc are approximately established by in- 
vestigations of the relic radiation (Spergel et al. 
2003) and the structure of the Universe at 2; <C 1 
detected in large redshift surveys (Percival et al. 
2001; Peacock et al. 2001; Efstathiou et al. 2001; 
Tegmark, Hamilton & Xu 2002; Verde et al. 2002) 
and weak lensing data (see, e.g., Hoekstra, Yee & 
Gladders 2002). The shape of the initial power 
spectrum on scales 10/i~^Mpc - l/i~^Mpc can be 
tested at high redshifts where it is not yet strongly 
distorted by nonlinear evolution (Croft et al. 1998, 
2002; Nusser & Haehnelt 2000; Gnedin & Hamil- 
ton 2002). Recent results on reconstruction of the 
initial power spectrum are summarized and dis- 
cussed in Tegmark and Zaldarriaga (2002), Wang 
et al. (2002) and Peiris et al. (2003). 

Analysis of the observed absorption spectra of 
high redshift quasars seems to be the most promis- 
ing method of recovering the power spectrum of 
initial perturbations at small scale. Indeed, the 
absorption lines trace the small scale hydrogen dis- 
tribution along the line of sight at redshifts z > 
2 when matter is not yet strongly clustered. The 
available Keck and VLT high resolution observa- 
tions of the Lyman— a forest provide a sufficiently 
rich database which can be analyzed with the help 
of statistical methods. 

The composition and spatial distribution of 
the observed absorbers is complicated and at low 
redshifts a significant number of stronger Ly-a 
lines and metal systems is associated with galax- 
ies (Bergeron et al. 1992; Lanzetta et al. 1995; 
Tytler 1995; Le Brune et al. 1996). However as 
was recently shown by Penton, Stock and ShuU 
(2002) , even at small redshifts some absorbers are 
associated with galaxy filaments while others are 



found within galaxy voids. 

These results suggest that the population of 
weaker absorbers dominating at higher redshifts 
can be associated both with the periphery of rare 
high density objects and with weaker DM pan- 
cakes formed by the non luminous baryonic and 
DM components in extended low density regions. 
They suggest that the Ly-a forest can be consid- 
ered as a low mass component of the generic LSS 
which is seen in simulated and observed spatial 
matter distribution. In turn, the relatively homo- 
geneous spatial distribution of absorbers implies 
a more homogeneous spatial distribution of both 
DM and baryonic matter as compared with the 
observed distribution of the luminous matter. 

Essential progress in understanding the process 
of formation and evolution of LSS has recently 
been achieved through numerous high resolution 
simulations which show that this process is mainly 
driven by the power spectrum of initial pertur- 
bations. These results were used in Demiahski 
& Doroshkevich (1999, 2002; hereafter DD99 & 
DD02) to consider the properties of the LSS in the 
context of the nonlinear theory of gravitational in- 
stability (Zel'dovich 1970; Shandarin & Zel'dovich 
1988). The statistical description of the process 
of structure formation and evolution developed in 
DD99 & DD02 links the directly measured charac- 
teristics of the structure such as the mass function 
of DM pancakes and walls, with the initial power 
spectum. 

This approach generalizes the Press-Schechter 
formalism (Press & Schechter 1974; Peacock & 
Heavens 1990; Bond et al. 1991; Loeb & Barkana 
2001) and describes the formation and evolution of 
all DM structure elements for the CDM or WDM 
initial power spectra without any smoothing or 
filtering procedures. The evolution of structure 
proceeds through random formation and merging 
of Zel'dovich pancakes, their transverse expansion 
and/or compression and subsequent transforma- 
tion into high density clouds and filaments. Later 
on the hierarchical merging of pancakes, filaments 
and clouds forms rich galaxy walls observed at 
small redshifts. The main stages of this evolution 
are driven by the initial power spectrum. 

As is well known, the Zel'dovich approach cor- 
rectly describes both the linear and mildly nonlin- 
ear stages of the structure formation but it cannot 
describe the processes of violent relaxation and the 
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final stage of formation of gravitationally bound 
object. In spite of this, the statistical approach 
proposed in DD99 and DD02 nicely describes the 
main properties of observed and simulated LSS 
(Demiahski et al. 2000; Doroshkevich, Tucker & 
Allam 2002). In particular, comparison of char- 
acteristics of simulated DM distribution and ob- 
served galaxy walls with theoretical expectations 
allows one to restore the amplitude and shape of 
initial power spectrum at appropriate scales. 

Early application of the truncated Zel'dovich 
approximation for the description of the Ly-a for- 
est was critically mentioned in Bond and Wads- 
ley (1997) but their criticism was not supported 
by more detailed investigation. In contrast, the 
Zel'dovich approach has been considered as very 
promissing by Hui, Gnedin and Zhang (1997) and 
Hui (1997) who applied it to describe the Doppler 
parameter and neutral hydrogen column density. 
However, these characteristics depend upon ran- 
dom variations of the gas entropy and, so, their 
theoretical description is quite problematic. 

The statistical approach proposed in DD02 was 
used in Demianski, Doroshkevich & Turchaninov 
(2003, hereafter Paper I) to describe properties 
of absorbers observed in 14 high resolution spec- 
tra. This approach deals with DM charactristics 
of absorbers. Evolution of DM component of ab- 
sorbers is driven mainly by the gravitational in- 
teraction and it can be approximately described 
theoretically. However in this approach the DM 
characteristics of absorbers are expressed through 
the observed ones what in turn requires a physical 
model of absorbers. 

Such model, introduced in Paper I, allows one 
to find a physically motivated combination of ob- 
served characteristics of absorbers which are con- 
sistent with theoretically expected ones for the 
DM pancakes. These results allows us to inter- 
pret this combination as actual characteristic of 
DM component of absorbers and confirm the self 
consistency of the proposed physical model of ab- 
sorbers. They show also that the statistical ap- 
proach proposed in DD99 and DD02 provides a 
reasonable description at least for so defined DM 
characteristics of absorbers. However, we cannot 
theoretically describe more complex evolution of 
the gaseous component of absorbers, their Doppler 
parameters and column density of neutral hydro- 
gen. 



Following Paper I, also in this paper we con- 
sider that absorbers are dominated by long-lived 
gravitationally bound and partially relaxed ob- 
jects formed in the course of both adiabatic and 
shock compression. Observations of galaxies and 
quasars at 2; > 3 and the reheating of the Universe 
at redshifts z > 6 demonstrate the importance 
of nonlinear processes at such redshifts and sup- 
port our approach. Numerous simulations (see, 
e.g. Frenk 2002) confirm the formation of high 
density DM filaments and sheets at similar red- 
shifts {z > 2). The typical size and mass of such 
objects progressively increase with time and now 
richer objects are seen as the LSS in the galaxy 
distribution. In this paper we assume that ab- 
sorbers at high z trace the DM structure which is 
qualitatively similar to the rescaled one observed 
at small redshifts. However, for less massive pan- 
cakes properties of the baryonic and DM compo- 
nents differ due to the infiuence of gaseous pres- 
sure. 

Here we use this approach to recover the initial 
power spectrum to unprecedently small scales. In 
contrast with previous investigations (Croft et al. 
1998, 2002; Nusser & Haehnelt 2000) we analyse 
the DM mass function of absorbers rather than 
the flux or smoothed density field. This means 
that our results are not restricted by the standard 
factors such as the Nyquist limit, the impact of 
nonlinear processes, the unknown matter distribu- 
tion between absorbers or their peculiar velocities. 
The reliability and precision of our method have 
been tested on the simulated DM and observed 
galaxy surveys (Demianski et al. 2000; Doroshke- 
vich, Tucker & Allam 2002). 

At scales > 1/i^^Mpc our estimates of the ini- 
tial power spectrum are consistent with the CDM- 
like one with precision ^ 10%. However, we found 
some evidence that at scales < 1/i^^Mpc the ini- 
tial power spectrum differs from the CDM-likc one 
suggesting a complex inflation with generation of 
excess power at small scales. Such excess power 
accelerates the process of galaxy formation at high 
redshifts and shifts the epoch of reheating of the 
Universe, which perhaps began at z = 20 ± 10 
(Kogut et al. 2003) and continued up to z w 6 
(Djorgovski et al 2001; Fan et al. 2001, 2003). 
More detailed discussion of the process of reion- 
ization can be found in Ciardi, Ferrara & White 
(2003) and Cen (2003). 
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At present we have only limited information 
about the properties of the background gas and 
UV radiation (Scott et al. 2000; Schaye et al. 
2000; McDonald & Miralda Escude 2001; Mc- 
Donald et al. 2000, 2001; Theuns et al. 2002 
a, b) and therefor some numerical factors in our 
model remain undetermined. This means that our 
approach should be tested on representative nu- 
merical simulations which more accurately follow 
the process of formation and disruption of pan- 
cakes and filaments and provide a unified picture 
of the process of absorbers formation and evo- 
lution (sec, e.g., Weinberg et al. 1998; Zhang 
et al. 1998; Dave et al. 1999; Theuns et al. 
1999). However, as was shown by Meiksin, Bryan 
& Machacek (2001), available simulations repro- 
duce quite well the characteristics of the flux but 
cannot restore other observed characteristics of 
the forest. This means that first of all numerical 
simulations should be improved (see Sec. 9.4). 

This paper is organized as follows. In Sees. 2 & 
3 a short theoretical description of structure for- 
mation and the physical model of absorbers are 
presented. In Sec. 4 the observational databases 
used in our analysis and redshift evolution of the 
observed absorber characteristics are described. 
Results of statistical analysis are given in Sees. 
5 and 6 and their reliability and precision are dis- 
cussed in Sec. 7. In Sec. 8 we link the mea- 
sured correlation function with the power spec- 
trum. Discussion and conclusions can be found in 
Sec. 9. 

2. Initial power spectrum and characteris- 
tics of DM pancakes 

2.1. Basic cosmological model 

In this paper we consider the spatially flat 
ACDM model of the Universe with the Hubble pa- 
rameter and mean matter density given by: 



We assume that the baryonic background is 
characterized by the following parameters 

{nb{z)) = 2.4 • 10-^cm-3(l -|- zf{nbh?lQm) , (2) 

Tig « 1.6 • lO^K, bbg = y^^^ « 16km/s . 

Here ilb, T^ghb^g are the dimensionless density, 
temperature and Doppler parameter of the bary- 
onic component, fc^ & ra}j are the Boltzmann con- 
stant and the mass of hydrogen atom. 

2.2. The initial power spectrum 

As a reference power spectrum of initial pertur- 
bations we take the standard CDM-like spectrum 
with the Harrison - Zel'dovich asymptotic, 

m = ^T\k/k,)Dw{k)Dj{k), fco = ^, (3) 

where A is the dimensionless amplitude of per- 
turbations, k is the comoving wave number. The 
transfer function, T{k), and the damping factor 
Dw{k) describing the free streaming of DM par- 
ticles were given in Bardeen et al. (1986). The 
factor Dj{k) describes the dumping of small scale 
baryonic perturbations (Jeans scale) and the bias 
between perturbations of DM and baryonic com- 
ponents (see, e.g., Mattarese & Mohayaee 2002). 
However, characteristics of the DM structure ele- 
ments are defined by (3) with Dj{k) = 1. 

For WDM particles the dimensionless damping 
scale, Rf, and the damping factor, Dw, are 

•'^ " 5 i Mdm ) ' ^ " exp[-r]Rf-{r]Rf) ] , 

where r] = k/ko and Mdm is the mass of WDM 

particles (Bardeen et al. 1986). The Jeans wave 
number, kj, and the damping factor, Dj, can be 
taken as 



H^{z)=H^nra{l + zf[l + n^/nmil + z)-^], (1) kj^^0.7bbg{l + z)H-\z), Dj^il + rjW)~\ 

iPmiz)) = ^n^{l+zf, Ho = 100/ikm/sMpc. rjj = ko/kj « 2 • IQ-^y^jyill^^^ . 

Here = 0.3 k^A = 0.7 are dimensionless mat- Here as a reference we use the values of the back- 

ter density and the cosmological constant (dark ground parameters specified by (2). 

energy), and h = 0.7 is the dimensionless Hubble ^j^^ amplitude of initial perturbations, A, is 

constant. measured reasonably well with different meth- 
ods. It is simply linked with erg, the variance of 



4 



the mass within a randomly placed sphere of ra- 
dius Sft-^Mpc, 

ai = 9A^ r d,,^T\n) ( ^^^^-r^^^A ' , (4) 
Jo \ Vs / 

where rj = k/k^ and 773 = S'q^lmh. Here we use 
the latest estimates (Spergel et al. 2003) for the 
model (1) 

0-8 « 0.055A « 0.9 ± 0.1, A « 16.4 ±1.82. (5) 

2.3. Coherent lengths and correlation 
functions of the initial perturbations 

For the spectrum (3) the coherent lengths of 
velocity and density fields, lyklp, are expressed 
through the spectral moments, m_2&mo, as fol- 
lows: 



ly 



m_2 



Too 



f 

Jo 



Mpc « 33.8ft-iMpc7?^ , (6) 



drir]T\ri)Dw{v)Dj{v) « 0.023, 



m 



drj7j'T'ir])Dw{v)Dj{v), Qo = 5- 
Too 



The analysis of the redshift distribution of 
Lyman-a clouds (Paper I) allowed us to estimate 
the moment Too that in turn restricts the expected 
Jeans scale and the mass of the dominant fraction 
of DM particles. These estimates 

go « (0.6 - 1.2) • 10-2, « 0.01 - 0.02 , (7) 

Mdm « (1.5 - 5)keV, Rj w (7.7 - 1.5) • 10"^ , 

are close to those of Narayanan ct al. (2000), 
Mdm > 0.75 keV, and Barkana, Haiman & Os- 
triker (2001), Mdm > 1 - 1-25 keV. However, as 
was discussed in Paper I, the Jeans damping, Dj, 
restricts abilities of this approach, and to test the 
small scale power spectrum we have to use charac- 
teristics of DM component of absorbers for which 
Dj = l. 

As was demonstrated in DD99 and DD02, 

main statistical characteristics of structure arc ex- 
pressed through the normalized longitudinal cor- 
relation function of the initial velocity field 



((q-v(qi))(q- v(q2))) 



q=^. (8) 



Here qi & q2 are the unperturbed coordinates of 
two particles a.t z = 0, q = |q|, and is the ve- 
locity variance. This function is expressed through 
the power spectrum: 



TO_2 



i dr^jf cosx -^T^y)Dw{y)= (9) 

Jt] y 




drj 



where x = grj, Ji/2{x) s-^d J5/2{^') ^rc Bcssel func- 
tions. From (9) it follows that ^^(0) = 1 and 
/o°° dq ^{q) = 0. For the CDM - like spectrum 
(3) and for the most interesting ranges go < IQ-^, 
0.5 > q > 10""^, the velocity correlation function 
can be fitted as follows: 



= CCDM 



1 - 



qo + 1.5uo{q) 



^/uo{q)ui{q) qo + Uo{9) 



(10) 



+ ql, ui = + 3 • IQ- 



Further on, we will use this function with go = 
0.01 as the reference one and will compare it with 
observational estimates of ^t,(g). This value of go 
lies in the range of (7) and corresponds to the mass 
of WDM particles Mwdm ~ 3 keV (or the Jeans 
damping scale for T^g w 10*K). 

2.4. Characteristics of DM pancakes 

We use the Zel'dovich approximation (Zel'dovich 
1970; Shandarin & Zel'dovich 1989) to relate the 
characteristics of the initial density perturbations 
with the properties of observed pancakes. As was 
discussed in DD99, DD02 and Paper I. three quan- 
titative characteristics of pancakes can be used for 
the comparison with observations. Two of them 
are the differential, A^^, and cumulative, Wq. dis- 
tribution functions (PDFs) of the dimensionless 
DM surface density (or column density) of pan- 
cakes, g. The DM surface density is defined as 
the dimensionless unperturbed distance at red- 
shift z = Q between DM particles bounding the 
pancake (8). The third characteristic is the di- 
mensionless mean number density of absorbers 
(Lyman-a clouds), riabsiQthr, z). All these char- 
acteristics depend upon the distribution function 
of the initial density perturbations which can be 
taken as Gaussian (Komatsu et al. 2003), upon 
the correlation function of the initial velocity field, 
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(.viq), and the survival probability of pancakes, 
Wslq), (DD99, DD02). The survival probability 
of clouds was introduced by Peacock fc Heavens 
(1990) and Bond et al. (1991) to describe, in 
the framework of Press-Schechter approach, the 
difference between the fraction of matter accumu- 
lated by clouds and their mass function. 

As was shown in DD99 and DD02, both PDFs 
and the survival probability depend upon the self 
similar variable 



c 



4r2(z)[l-e„('?)] 



(11) 



where the 'time' t{z) describes the growth of per- 
turbations due to the gravitational instability. For 
the ACDM cosmological model (1), and for ^; > 2 

we have 



t{z) w To 



1 + 1.2nr. 

2.20777, 



1/3 



1.27T0 

1 + z ^ 1 + z 



To 



= Am_2/\/3 » 0.22 ± 0.02 , 



, (12) 



(13) 



where A was introduced in (3) and (5). 

For the Gaussian initial perturbations the mat- 
ter fraction accumulated by structure elements is 
proportional to erfc((^)=l-erf((^), where erf((^) is 
the standard error function. Following Peacock & 
Heavens (1990) and Bond et al. (1991), we as- 
sume that the mass function of the DM column 
density of pancakes and the matter fraction accu- 
mulated by pancakes differ by the survival prob- 
ability, Ws{C)- For low mass pancakes Wg <C 1, 
what reflects stronger correlation of small scale 
perturbations, and Wg — > 1 with increasing mass 
of pancakes. 

Analysis of simulated and observed structure 
at small redshifts (DD99; Doroshkevich, Tucker & 
Allam 2002) shows that for richer pancakes and 
galaxy walls the cumulative and differential PDFs, 
Wq & Ng, for the DM column density can be writ- 
ten as follows: 

with the survival probability Wg = erf(-\/C). How- 
ever, other expressions for the survival probability 
can be also used and, for comparison, here we con- 
sider the model with Wg = erf^(VC) and 



This problem is further discussed in Sec. 7. 

The dimcnsionless mean number density of ab- 
sorbers, nabs{qthr, z), wlth the above assumptions 
is given by: 

c erfc(Cthr)(l + ^)^ 



Holy {q[Cthr)) 

Here erfc{(thr) describes the progressive matter 
concentration within pancakes with ( > (thr = 
C{z, qthr) and quir is the minimal (threshold) DM 
column density of pancakes in the sample under 
consideration (DD02, Paper I). For Ws =e^:i{^/C,) 
and for qthr » go, we get 



{qiCthr)) = 



1 + 



4\/7rCt/irerf(yCt^) + 2 exp(-Ct/ir)" 



7rexp(Ct;ir)[l - erf^(yCt^)] 

where C,thr{z) = qthr I ^t'^ [z) . Similar relations can 
be also written for other survival probabilities. 

For pancakes with small qthr ^ ?o the red- 
shift variations of the mean number density 
nabs{z,qthr) arc described by the expression 



cV3 (1 + 2)2 exp(-C 



thr) 



WsiUr), (17) 



with go as given by (7) (DD02, Paper I). This 
approach considers only gaseous pancakes and, 
as was noted in Sec. 2.3, its sensitivity to the 
small scale perturbations is restricted by the Jeans 
damping. 

2.5. Reconstruction of the initial power 
spectrum 

Using the measured rcdshift, z, and DM column 
density of absorbers, q, we determine the cumu- 
lative PDF of absorbers Wo6s[> q/T^z)] and, for 
each q/T{z), we compute (17) and a"^ = (q^) — (q)"^. 
For a chosen Wq{C,), we solve numerically the equa- 
tion 



Wobs[> q/T{z)] = Wq{C), 



(18) 



with respect to C{q,T) and, thus, we obtain the 
function 

„2 



Uq) = 1 



(19) 



dq 



4r2C ■ 

The observed functions Nq{Q and nabs{z) can be 
compared with corresponding expectations (14) - 
(16) for the correlation function £,v{q) (19). This 
S,v{q) can be also compared with the reference 
function ^cDM{q) (10). 
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3. Determination of the DM column den- 
sity 

To connect the main observational characteris- 
tics of the Lyman-a absorption Unes, namely, the 
redshift, Zabs, the column density of neutral hy- 
drogen, Nhi, and the Doppler parameter, h, with 
the dimensionless DM column density, q, wc use 
the improved model of absorbers discussed in Pa- 
per I. Here we only briefly summarize the main 
assumptions and relations. Detailed discussion of 
the model can be found in Paper I. 

3.1. Physical model of absorbers. 

Following Paper I we assume that: 

1. The DM distribution forms an intercon- 
nected structure of sheets (Zel'dovich pan- 
cakes) and filaments. The baryonic compo- 
nent follows the DM distribution and there- 
fore the main parameters of absorbers are 
approximately described by the Zel'dovich 
theory of gravitational instability (Sec. 2; 
DD99; DD02). 

2. The majority of DM pancakes are partly re- 
laxed, and long-lived, but the observed prop- 
erties of absorbers vary with time due to the 
successive merging, transversal compression 
and/or expansion and disruption of DM pan- 
cakes. 

3. Gas is trapped in the gravitational potential 
wells formed by the DM distribution. The 
gas temperature and the observed Doppler 
parameter, b, trace the depth of the DM 
potential wells. We consider the possible 
macroscopic motions within pancakes as 
subsonic and assume that they cannot es- 
sentially distort the measured Doppler pa- 
rameter. 

4. The gas is ionized by the UV background 

and for the majority of absorbers ionization 
equilibrium is assumed. 

5. For a given temperature, the gas density 
within the potential wells is determined by 
the gas entropy created during the previ- 
ous evolution. The gas entropy is changing, 
mainly, due to shock heating in the course 
of pancakes merging, bulk heating by the 



UV background and local sources, and due 
to radiative cooling. Random variations of 
the intensity and spectrum of the UV back- 
ground enhance random scatter of the ob- 
served properties of absorbers. 

This model implies that the hy-a clouds trace 
the LSS in the DM component which is qualita- 
tively similar to the rescaled structure observed at 
2; < 1 in the galaxy distribution. It links the char- 
acteristics of the DM structure and the observed 
absorption lines and, as was shown in Paper I, pro- 
vides a reasonable self consistent description, at 
least, of richer absorbers. Our model considers dis- 
crete absorbers and, so, it essentially differs from 
the popular model of " fluctuating Gunn-Peterson 
approximation" discussed more thoroughly in Sec. 
9.4. 

Our approach cannot be applied to objects for 
which the gravitational potential and the gas tem- 
perature along the line of sight strongly depends 
on the matter distribution across this line. How- 
ever, the anisotropic halos of fllaments and clouds 
with moderate Nhi can be considered as 'pancake- 
like'. 

Special attention require description of weaker 
absorbers with h < b^g. This subpopulation con- 
tains "artificial" caustics (McGill 1990) and ab- 
sorbers identified with slowly expanding under- 
dense regions (Bi & Davidsen 1997; Zhang et al. 
1998; Dave et al. 1999). Such absorbers produce a 
short lived noise, which is stronger at higher red- 
shifts 2; > 3. However, such absorbers can also 
be formed within low temperature regions arising 
owing to spatial nonhomogeneities of the UV back- 
ground and expansion rate. Strictly speaking our 
model does not include such absorbers because we 
assume that b > b^g for all absorbers formed due 
to compression of background matter and we sup- 
pose that bi)g is constant and it does not depend 
on the position or time. However, we can include 
absorbers with b < btg into consideration assum- 
ing that they arc adiabatically compressed but not 
bound gravitationally. 

3.2. DM column density of absorbers 

In this section we use the physical model dis- 
cussed in Sec. 3.1 to introduce the basic relations 
between theoretical and observed characteristics 
of pancakes. The most important characteristics 
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of DM pancakes are presented (without proofs) as 
a basis for further analysis. More details are given 
in DD99, DD02, and Paper I. 

The fundamental characteristic of DM pancakes 
is the dimensional, i^{z), or the dimensionless, q, 
DM surface density (the dimensionless DM col- 
umn density) : 



{Pm{z))lvq 



(20) 



where the coherent length of the initial velocity 
field ly was introduced in (6). As was discussed in 
Sec. 2.4, lyq is defined as the unperturbed distance 
at rcdshift = between DM particles bounding 
the pancake. 

As was assumed in Sec. 3.1, we consider ma- 
jority of absorbers as gravitationally bound and 
relaxed long lived DM pancakes. Their properties 
are slowly changing mainly owing to the process 
of merging and transverse expansion and compres- 
sion. For such absorbers, the gas temperature. 



T„, and the Doppler parameter, b cx 



(2) are 



closely linked with the gravitational potential of 
the DM pancake and we get: 



P - u2 ~ i + _ - 
% 



5 blM^))6r, 



l + 5o|-(l + 0),(21) 



So=-: 



3 / HqI^ 



0.13 /16km/s 



\ bbgiz) 



10 V hi 

Here Sm is the mean overdcnsity of the DM pan- 
cake and the random factor 0$ ^ 1 describes 
the nonhomogeneity of DM distribution across the 
pancake. 

Under the condition of ionization equilibrium, 
for the column density of neutral hydrogen, we get 

Nhi = ^^^^^hq(3-'/'{l + zf (22) 

_2-10^^cm-^ q^{l + zf f3^ Sb 
Gl2{z) /37/2 /32_i^' 



ri2(-z) 008 61 / Qmh^ 0.02 y / b, 



Ogiz) 0.5 V 0.13 fib/iV VlSkm/s 



Obg 



7/2 



(23) 



where = T^/IO" 



F-Y characterizes the 



rate of ionization of hydrogen by the UV back- 
ground, Sb is the overdensity of baryonic compo- 
nent, cos^ describes a random orientation of the 



absorber and the line of sight and the function 
Qg{z) ~ 1 describes the action of several random 
and unknown factors (more details are given in 
Paper I). 

The equations (21) and (22) relate three inde- 
pendent variables, namely, q.SbSzSm- To find the 
DM column density, q. it is therefore necessary to 
use an additional relation which connects the ba- 
sic properties of absorbers. In Paper I we assumed 
that (3 > 1, Sm Sb- These assumptions are valid 
for richer absorbers formed due to strong shock 
compression when the essential growth of the en- 
tropy and temperature suppresses differences be- 
tween the DM and baryonic components and so 
their distribution across the absorber is similar. 
For the DM column density of such absorbers we 
get 



N, 



HI 



2 • 10i6cm-2 



^3/2(^2 _ ^) 
(1 + Z)6 



Gi2{z). (24) 



This relation describes correctly properties of hot 
absorbers but decreases the column density of rare 
absorbers with 5m < Sb- 

However, formation of weaker DM absorbers 
with /? ~ 1 is accompanied by adiabatic infiow 
of baryons into the potential well. In this case, in 
accordance with the polytropic index 'jb = 5/3 for 
baryons, we can expect that 



9$. 



So 



/32- 1 



>Sb, 



(25) 



> 



/33/2(/32_ 1)1/2^-1/2(1 _^^)-l/2 



where Sq was introduced in (21). This expression 
quantifies the bias between the compressed DM 
and baryonic components arising due to the im- 
pact of gaseous pressure (Jeans damping). For 
such absorbers we get from (22): 



HI 



G12S0 



2- 10l6cm-2 J /J3/2(l-K0)5 ' 



(26) 



Both expressions, (24) and (26), define the di- 
mensionless column density of DM component cor- 
rected for the impact of gaseous pressure. These 
relations can be successfully applied to absorbers 
formed due to adiabatic and strong shock com- 
pressions. However, the boundary between these 
limiting cases must be established a priori- So, to 
discriminate absorbers described by (24) and (26), 
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we can use a threshold Doppler parameter, Pthr, 
which, in fact, characterizes the Mach number of 
the inflawing matter. Thus, the weak adiabat- 
ically compressed absorbers can be conveniently 
separated with a threshold, 

l</3</?t/.r«1.5-2, (27) 

under condition (25) which can be reformulated 
as a restriction on the column density of neutral 
hydrogen, Nhi- Absorbers with /? > Pthr are con- 
sidered as dominated by shock compressed gas and 
their DM column density is given by (24). 

As was noted above, we cannot correctly de- 
termine properties of absorbers with h < bbg in 
the framework of our model because they violate 
the relation (21). However, we can consider them 
as adiabatically compressed, short lived unrelaxed 
and gravitationally unbounded pancakes. Under 
these conditions absorbers with b < bbg can be 
described by the expression (26) which does not 
depend on bbg. 

As was found in Paper I, we can approximate 
the combined action of redshift variations of the 
UV background and random factors in expressions 
(24 & 26) by a two parameters relation 

{Gi2{z)) ^ Go{l + a^z) , (28) 

where Go and are constants. The precision of 
this approach is moderate but a suitable choice 
of Go and provides a reasonable statistical de- 
scription for majority of absorbers. 

3.3. Statistical characteristics of absorbers 

The model of absorbers discussed in this section 
is in fact independent from the statistical charac- 
teristics discussed in Sec. 2 . It introduces the 
physically motivated combination of observed ab- 
sorber characteristics, q, as given by (24) and (26), 
which can be compared with theoretical expecta- 
tions introduced in Sec. 2 . As it is shown in 
Sec. 5, the fimctions (il^ z) (11) and Wg found 
with this approach are consistent with expecta- 
tions what verifies the interpretation of q as the 
DM column density of absorbers. 



4. Mean characteristics of the observed 

sample 

4.1. The database 

In our analysis we use 14 high resolution QSO 
spectra listed in Table I . To decrease the scatter, 
we selected the sample of absorbers requiring that 
IQi^cni-^ < Nhi < lO^-'^cm-^ & < 90 km/s. 

Absorbers with b < bbg must be excluded from 
the sample because for them the overdensity Sb = 
(b/bbg)"^^^ < 1 . For bbg = 16km/s our sample in- 
cludes 3752 absorbers at redshifts 1.7 < z < 4. 
Subpopulation of 500 absorbers with b < bbg 
contains "artificial" caustics (McGill 1990) and 
absorbers which are most probably situated within 
regions where the average density and tempera- 
ture are smaller than the values specified in (2) 
(Bi & Davidsen 1997; Zhang et al. 1998; Dave 
et al. 1999). As was noted above we can include 
them into consideration assuming that they are 
short lived gravitationally unbound pancakes. 

For all models discussed in this paper we de- 
scribe the redshift variations of the UV back- 
ground intensity by the relation (28) with 

Go = 12 - 15, a-y = -0.22 . (29) 

Estimates of the intensity of UV background 
through the proximity effect (Scott et al. 2000) 
give «2±lat2~3. Similar estimates were 
obtained in McDonald & Miralda-Escude (2001). 
These results are consistent with (23), (28) and 
(29) when 

(Qg) ~ 0.4(66g/16km/s)-^/2 . 

Such a choice is in accordance with expectations of 
our model of absorbers and results in weak redshift 
variations of the function {({z)). 

4.2. Redshift variations of absorbers char- 
acteristics 

Here we consider the redshift variations of the 
observed parameters of absorbers, namely, the 
neutral hydrogen column density, logNni, the 
Doppler parameter, b, and distances between ab- 
sorbers 

= ^ - 5.5 . ^O^^^^'^-'Mpo.iSO) 
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For our sample of absorbers, the redshift varia- 
tions of the mean observed parameters and the 
mean DM column density measured by the self 
similar variable, {C{z)) (11), are plotted in Fig. 1 
together with approximate fits : 

{C{z)) w 0.73 ± 0.08, {b{z)) w (31.9 ± 1.6)km/s , (31) 

{logNHi) « 13.2±0.22, {Dsep) « {2.5-0.6z)h-^Mpc. 

These results are consistent with the main results 
obtained in Paper I. It is interesting to notice that 
the redshift variations of {({z)) and {b{z)) plotted 
in Fig. 1 are clearly correlated. 

The PDFs Pb, Phi and Psep, for the Doppler 
parameter, b, the hydrogen column density, Nhi, 
and the absorbers separation, Ds^p, are plotted in 
Fig. 2 for three redshift intervals, namely, 1.7 - 
2.5 (1357 lines), 2.5 - 3 (1235 lines) and 3 - 4.5 
(1160 lines). This Fig. indicates that the shapes of 
PDFs weakly vary with redshift, what is consistent 
with the assumption that absorbers are long-lived 
and relaxed. 

The mean values of the measured (b) and {Nhi) 
for these intervals of z are quite stable: 

(6) = 28.3km/s, 31.1km/s, 30.0km/s, 

(Nhi) = lO^^-^cm-^ 1013.6^^^-2^ 10^3.8^^-2 ^ ^32) 

(£>^ep) = 1-75/i-iMpc, 1.34/i-iMpc, O.Mh-^Mpc. 

To quantify the redshift variations of the PDFs 
we calculate the mean values and dispersions of 
the relative differences, 

c^ = [Pl'^-Pt]/[Pl^+Pt], 

between PDFs for the full sample (-P/*"* , 3752 lines) 
and for each individual redshift interval {P^, k=l, 
2, 3), where i = 1,50 numbers the bins with a 
step of 0.1 of the mean value. For the most repre- 
sentative first 25 bins (cf ) = with the standard 
deviation ctc ^ 0.1 — 0.15. 

These results indicate also the homogeneity of 
our sample. They confirm that the redshift evo- 
lution of the observed parameters of absorbers is 
very similar to the evolution of their mean val- 
ues (31). It is an evidence in favor of the model 
of long-lived gravitationally confined absorbers. 
More detailed discussion of the observed spectra 



can be found in Kim, Cristiani & D'Odorico (2002) 
and Kim et al. (2002). 

It is interesting to estimate also the mean 
matter fraction accumulated by absorbers, (/m). 
Rough estimates of fm can be found by compar- 
ing the unperturbed mean size, lv{q), with the 
mean separation of absorbers, D/Nabs for each 
individual spectrum, we get 

ifm) ~ Nabslv{q)/D « 0.87 ± 0.14 . (33) 

Here Nabs and D are the number of absorbers in 
the spectrum and its length along the line of sight. 
Averaging is performed over all fourteen spectra in 
our data base. 

This value (fm) exceeds the theoretical expec- 
tation (fm) ~ 0.5 — 0.6 at 2; ~ 3 (Paper I) what 
indicates a limited precision of this approach. In 
spite of this, it verifies that, in accordance with 
the physical model, at redshifts z < 3 at least 
half of matter could be accumulated by absorbers. 
This estimate does not contradict the observed 
Gunn-Peterson effect in Hell at 2: ~ 3 (Jacobsen 
et al. 1994) what indicates that essential fraction 
of baryons is homogeneously distributed. However 
quantitative estimates of this fraction strongly de- 
pends upon the intensity of UV background ioniz- 
ing Hell. 

5. Correlation function of the initial veloc- 
ity field 

In this Section we consider the cumulative 
and differential distribution functions, Wq{Q and 
Nq{(^), and the correlation function of the initial 
velocity field, which can be found from Wq{(^) 
as was described in Sec. 2.5. These functions 
depend upon several parameters which are known 
with limited precision. To test their influence we 
consider below three models of absorbers. 

5.1. Three models of absorbers 

Our results weakly depend upon the parameter 
Pthr introduced in Sec. 3.2 to separate absorbers 
formed by adiabatic and shock compression. The 
choice of jSihr is restricted by the shape of the fimc- 
tion S^y at larger q (10) as Pthr < 2. In our model 
the amplitude of initial perturbations and 

To = 0.22, Ahr = l-5, bbg = mkm/s, (34) 
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remain fixed. As is seen from (32), Pthr ^ {b)/bbg 
and it turns out that ~ 700 — 1 ODD of absorbers 
in our sample have been adiabatically compressed. 
They form the left side of the PDF Pt plotted in 
Fig. 2 while richer absorbers with /3 > Pthr ~ 
(b) /bbg have been presumably formed by merging 
and shock compression. 

The cumulative and differential distribution 
functions depend also upon the survival proba- 
bility, W, =erP(yC), (14, 15) described by the 
power index n. Finally, of course our results de- 
pend upon the number of absorbers in the sample, 
Nsmp, used in the analysis. In particular, in spite 
of the relatively small number of absorbers with 
b < bbg when they are included they noticeably 
change the correlation function 

To test the influence of these factors we consider 
three models with parameters: 

n = l, iV,„p = 3 752, Go = 12.5, (35) 
n=l, iV,„,p = 4 248, Go = 12.5, (36) 

n = 2, 7V,„p = 4 248, Go = 15.5, (37) 

for the Model 1, Model 2 and Model 3, respec- 
tively. Parameter Go was introduced in (28) to 
characterize the intensity of the UV background. 
The Model 1 includes only absorbers with b > bbg 
whereas Models 2 & 3 include also absorbers with 
b < bbg and differ by the choice of survival prob- 
ability, Ws, and the intensity of the UV back- 
ground. 

5.2. Cumulative PDFs for dimensionless 
DM column density 

As was discussed in Sec. 2.5, the observed cu- 
mulative distribution function, Wq{> q/r), can be 
estimated using q/r obtained from our sample of 
observed absorbers with the model described in 
Sec. 3 and parameters (34-37). For these mod- 
els the correlation function is well fitted by the 
expression 



^fit = 1 - 1.5gJ 4x4' 0.5 > g > 10-3 , (38) 

Pl = 1.3•10-^ J^l = 10-^ J^l=2•10-^ 
P2 = 7 • 10-^, P2 = 4 • 10-^ P2 = 2 • 10-3 , 



for the Model 1, Model 2, and Model 3, respec- 
tively. The function (38) becomes identical with 
the reference correlation function (10) ^cdm for 
q>Pi> P2- 

The measured functions Wq for the three mod- 
els are plotted in Fig. 3 together with fits (14) 
and (15) for given by (38) and with a scatter 

2/y/Nabs where Nabs is the number of absorbers 
in each bin. The diflference between the measured 
and fitted functions Wq is characterized by the re- 
lations: 

r (0.96±G.2)erf2(yC), = 0.5 , 
Wq^l (1.0±0.1)erf2(yC), X'=0.2, (39) 
[ (l.l±0.2)erf3(^), x' = 0.5, 

where ( is given by (11) and (38). is found for 
ratios of measured and fitted Wq/erP+^(-\/C) for 
26 measured points. For ^„ given by (38) all fits 
are reasonably consistent with the measured Wq. 

5.3. Differential PDFs for dimensionless 
DM column density 

For all three models the differential PDFs were 
found from measured q and r with given by 
(38). These PDFs are quite sensitive to the func- 
tions and Ws, what allows one to check the 
correctness of the choice of these functions. 

The differential PDF, Nq{C) = {C)dWq/dC, is 
plotted in Fig. 4 for the Model 1, Model 2, and 
Model 3 together with the fit 

„.,0»(!i±S<<>exp,-0^. (40) 

This expression fits well the observed PDF for 

(W« (0 = 0.85, 0.87, 1.15, (41) 

for Models 1, 2, & 3, respectively. Here {(th) = 
0.82 & 1.05 are the mean values for A^^ (40) with 
n = 1 and n = 2. For ( < 0.2(^) the measured 
PDF becomes sensitive to the possible incomplete- 
ness of the sample. 

The concordance of the measured and expected 
PDFs can be also checked by the Kolmogorov- 
Smirnov test. For the three models plotted in 
Fig. 4 the cmnulative PDFs are consistent with 
erf'+^(v^) with probabilities 

Pks = 0-7S, Pks = 0-76, Pks=0A7, (42) 
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respectively. These results confirm, for a given set 
of DM column density, a close link between the ac- 
cepted survival probability, Wg , and the measured 
correlation function, ^^j. For ( > {() the expo- 
nential shape of the function Nq coincides with its 
theoretically expected shape for Gaussian initial 
perturbations and the CDM-like power spectrum 
(10), (11), (14). 

5.4. Correlation function of the initial ve- 
locity field 

As was described in Sec. 2.5 we can derive 
the correlation function of the initial velocity field, 

£,v{q), by solving Eq. (18). This function depends 
on both the survival probability, Wg, and the sam- 
ple of absorbers (see 34-37). 

For all three models, the measured correlation 
functions (q) are plotted in Fig. 5 together with 
£,fit (38) and with the correlation function ^cdm 
(10). For all models variations of the ratios of 
measured and fitted functions, (1 — iv)/{^ — ^fit), 
are quite moderate and are quantified as follows: 

(1 - e.)/(l - ~ 1 ± 0-2, 0.3 > g > 10-3 , (43) 

X^ = 0.9, 0.4, 0.7, 

for the Model 1, Model 2 and Model 3, respec- 
tively. Here are found from 26 measured points. 
This fact confirms the self consistency of our mod- 
els and the choice of parameters (29) and (34-37). 

At the same time, at larger scales, 0-3 > q > 
0.03, 10/i-iMpc> lyq > l/i-^Mpc, for all three 
Models the differences between the measured and 
CDM- like correlation functions are negligible 

(1 - C,)/(l - ^cdm) « 1 ± 0.15 , x' = 0.2 , (44) 

where are found for the ratios (1 — Cv)/{^ — 
^CDm) for 14 measured points. These results 
demonstrate that for such scales the differences 
between the CDM -like initial power spectrum and 
the measured one do not exceed 15%. 

However, at small scales all measured correla- 
tion functions differ from the ^cdm (10). Quanti- 
tatively these differences are characterized by the 
XcDM found for the ratios (1 - C«)/(l - £,cdm) 
for 26 measured points. For the Model 1, Model 2 
and Model 3 we get, respectively : 

XcDM = 550, Xcz?M = 120, Xcz3M = 33.(45) 



These results show that this difference between 
the measured and £,cdm decreases when for- 
mation of pancakes with small mass is suppressed 
more strongly i.e. when the survival probability is 
Ws =erf2(VC). 

6. Redshift distribution of absorbers 

The redshift distribution of absorbers only 
weakly depends on the properties of weaker ab- 
sorbers because they are relatively rare, only 
about ~ 20 % of absorbers have q < 0.03. It 
depends mainly upon measured redshifts of ab- 
sorbers and is only weakly sensitive to the physi- 
cal model of absorbers discussed in Sec. 3. This 
means that it is an important independent test 
of the model of absorbers formation considered in 
Sec. 2. 

The redshift distribution of absorbers is quite 
well fitted by the relation (16) which is applied to 
all samples of absorbers with qthr > qo- To illus- 
trate this property we plot in Fig. 6 the measured 
function {nabs{z)) for all 4 250 absorbers together 
with two fits (16). 

As was described in Paper I, for each redshift in- 
terval, Az = 0.1, all absorbers were organized into 
an "equivalent single field" by arranging them one 
after the other along the line of sight. The mean 
number density of absorbers, riabsj and its scat- 
ter were found by comparing the equivalent single 
field with a corresponding Poissonian distribution. 
The measured function Uabs is plotted in Fig. 6 
together with the best fits (16) for Model 1 and 
Model 3. Both fits are very similar to each other 
what demonstrates that the observed redshift dis- 
tribution of absorbers only weakly depends on the 
shape of the initial power spectrum. However, the 
observed redshift distribution of absorbers verifies 
the Gaussianity of the initial perturbations since 
it is manifested in the characteristics of the Ly-a 
forest. 

7. Reliability and precision of our esti- 
mates 

To measure the initial power spectrum it is nec- 
essary to use a model of absorbers which links the 
observed Ly-a clouds with the DM spatial distri- 
bution. The model introduced in Sees. 2 & 3 pro- 
vides a quite reasonable and self consistent sta- 
tistical description of the basic properties of ab- 
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sorbcrs (Paper I and Sees. 4&5). For the three 
models under consideration the measured corre- 
lation function of the initial velocity field, ^v, is 
consistent with the CDM like one for q > 0.03, 
what verifies the reliability of our approach. 

However, at small scales {q < 0.03) the mea- 
sured function differs from the CDM-like one. 
These differences are caused by the deficit of ab- 
sorbers with small g/r & g and, therefore, it can be 
artificially enhanced by an action of several fac- 
tors. The most important ones are the possible 
incompleteness of the observed sample of weaker 
absorbers, overestimates of the DM column den- 
sity by our model of absorbers (24, 26), and the 
choice of the survival probability (14, 15). 

Usual estimates of the possible incompleteness 
of the sample (< 20%) are based on the quality 
of the observed spectra but for compiled samples 
of absorbers such estimates are not reliable. The 
best way to test the influence of this factor is to re- 
peat the analysis with a richer sample of observed 
absorbers. 

By comparing the Model 1 with the Model 2 
we test the sensitivity of our results to variations 
of the observed sample. In spite of the fact that 
the condition b < btg is satisfied only by a rela- 
tively small fraction of absorbers (~ 500) their in- 
fluence is quite significant because when they are 
included they almost double the number of weaker 
absorbers. In particular, the difference between 
the measured and the reference one (10) is re- 
duced by half but it is not eliminated. As is seen 
from (45), this difference remains clearly notice- 
able even for the Model 3 for which both factors, 
larger number of absorbers and larger power index, 
n = 2, of the survival probability work together to 
decrease this difference. 

On the other hand, extension of the sample by 
adding 355 absorbers with lO^^cm"^ < Nhi < 
lO^^cm"^ from published spectra of four quasars 
(2126-158, Giallongo ct al. 1993; 1700-1-642, Ro- 
driges et al. 1995; 1225-1-317, Khare et al. 1997; 
1331-M70, Kulkarni et al. 1996) does not change 
our previous results (43, 45). This test demon- 
strates that the measured velocity correlation 
function strongly depends on the representativ- 
ity of tlu; sample of weak absorbers and estimates 
of their Doppler parameter but it only weakly 
depends on the sample of stronger absorbers. 



Our choice of the survival probability in the 
expression (14) is based on the measured distri- 
bution of walls in observed galaxy surveys and in 
simulations rather than on theoretical argmncnts. 
Comparison of results obtained with the expres- 
sions (14&15) illustrates the importance of this 
factor and shows that the measured velocity corre- 
lation function could be noticeably different from 
that derived in a standard way from the CDM- 
like power spectrum only for a moderate survival 
probability, Ws{q,C)- The analysis repeated for a 
wider set of functions Wg {q, () confirms that the 
measured functions Wq and Nq strongly restrict 
the acceptable shape of the survival probability, 
but for a suitable choice of Wg , the difference be- 
tween £,cDM and the measured velocity correlation 
function becomes moderate. 

These results show that the merging of low 
mass pancakes is an important factor and it should 
be thoroughly investigated. As the first step of 
such investigation the PDFs Nq were measured 
for three sets of DM walls selected at z = with 
a small threshold richness from two standard high 
resolution simulations (Cole et al. 1998; Jenkins 
et al. 1998). These PDFs are plotted in Fig. 7 
together with fits (14). As is seen from this Fig- 
ure, in all the cases there is a noticeable excess 
of low mass objects. It can be related to the ex- 
pected disruption of compressed DM walls into a 
system of high density clouds, what decreases the 
observed pancake column density. These results 
favor a moderate survival probability of low mass 
absorbers. However, such analysis should be con- 
tinued with more precise simulations. 

As was noted in White & Croft (2003) the for- 
mation of LSS elements is accompanied by their 
partial disruption and creation of low mass clouds. 
This process is more important for formation of 
galaxies and their low mass satellites what was 
also discussed in DD02 and Paper I. This process 
can increase the fraction of low mass absorbers in 
comparison with theoretical expectations, what, 
however, can only partially suppress distortions of 
the initial power spectrum considered in this pa- 
per. 

This discussion confirms that quantitative es- 
timates of the shape of the function ^„ at small 
scales actually depend upon the sample used in the 
analysis and the model linking the observed ab- 
sorber characteristics with their DM column den- 
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sity and the correlation function ^y. However, dif- 
ferences between the measured and the CDM-hke 
correlation functions remain quite significant for a 
wide range of models. 

By analyzing the distribution of Lyman-a 
clouds we can test the initial power spectrum of 
perturbations to unprecedently small scale and 
both outcomes confirmation of the CDM-like 
shape of the power spectrum or detection of dis- 
tortions are equally important. Results of such 
analysis can be improved with a more detailed 
numerical simulations and a richer sample of ob- 
served absorbers. 

8. The initial power spectrum 

For larger scales, q > 0.1, our approach de- 
scribes reasonably well the properties of richer 
absorbers and our analysis verifies the CDM-like 
shape of the correlation function ^„ and the initial 
power spectrum (3). The correlation function 
is well fitted by the one parameter relation (10), 
where the parameter go is expressed through the 
spectral moments, mo and rn_2 (6), and it de- 
pends upon the mass of dominant DM particles. 
For Mdm « 3 - 5 keV, go « (1 - 0.7) • 10"^ is 
expected. 

However, at small scales the shape of the ob- 
served correlation function differs from the 
CDM-like one (10), what indicates a more compli- 
cated shape of the initial power spectrum at small 
scale. The initial power spectrum can be restored 
from the function ^y. Using (9) we obtain 



P{r])cx [ dqq^iq) 
Jo 



sm{qr]) + 2 



cos{qr]) 



qt] 



(46) 



The precision and reliability of our estimates are 
however limited (Sec. 7) and we measure the func- 
tion ^(g) only in limited range of g. So, here we 
will present only examples of spectra which gen- 
erate the correlation functions ^t, similar to the 
measured one. As was noted in Sec. 3, the expres- 
sions (24) and (26) define the DM column density 
corrected for the .Jeans damping and therefore as 
a reference we consider the standard normalized 
spectrum (3) for the mass of dominant fraction of 
DM particles Mdm = 3 keV, go = 0.01: 

Q{r,)=rjT^irj)Dw{v), V = k/ko, Dj{r,) = 1 . {A7) 



Distortions of the spectrum (47) are described by: 
AQiv) = exp[-0.5(77 - 77p)V^p] , (48) 



Art 



271 dp 

5 • 10"- 



Ap = 2.8 • 10" 



rjp = 170, (Tp = 65 , 

rjp = 170, ap = 280 , 

for Model 1 and Model 3, respectively. The corre- 
lation functions S,v generated by the spectra (47) 
and (48) and these spectra are plotted in Fig. 8 
and Fig. 9 . 

The interpretation of these distortions is not 
unique because of very limited available informa- 
tion. We can estimate the function ^„ only at g > 
10~^ where its dependence upon the mass of DM 
particles is moderate and even for Mdm > 10'' 
keV it increases of about 1.5 times only. This 
means, for example, that even so essential increase 
of Mdm cannot provide noticeable effects and the 
most probable source of the observed distortions 
are adiabatic or, perhaps, isocurvature perturba- 
tions predicted by some models of the early Uni- 
verse and inflation. 

Recent WMAP measurements confirmed that 
as expected adiabatic perturbations dominate on 
large scale (Peiris et al. 2003) and that they are 
Gaussian (Komatsu et al. 2003). However, these 
results do not preclude the existence of distor- 
tions of the standard CDM-like power spectrum 
at small scales due to both adiabatic and isocur- 
vature perturbations (Peiris et al. 2003). 

The measured distortions of the power spec- 
trum can be considered as an observational ev- 
idence in favor of the one field infiation with a 
complicated infiation potential (see, e.g., Ivanov, 
Naselsky & Novikov 1994) or multiple fields in- 
flation (sec, e.g., Polarski & Starobinsky 1995; 
Turok 1996). Both models generate adiabatic or 
isocurvature deviations from the simple CDM-like 
power spectrum. More detailed discussion of such 
models can be found, for example, in Peiris et al. 
(2003). 

9. Summary and Discussion. 

In this paper we continue the analysis of a 
sample of ~ 4 500 observed absorbers initiated in 
Paper I and based on the statistical description 
of Zel'dovich pancakes (DD99, DD02). This ap- 
proach allows one to connect the observed charac- 
teristics of absorbers with fundamental properties 
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of the initial perturbations and demonstrates the 
common origin of absorbers and the Large Scale 
Structure observed in the spatial galaxy distribu- 
tion at small redshifts. 

The physical model of absorbers proposed in 
Paper I and essentially improved in Sec. 3 links 
three observed characteristics of the Lyman-a 
clouds, namely, the redshift, z, the Doppler pa- 
rameter, h, and the column density of neutral hy- 
drogen, Nhi, with the DM and baryonic column 
densities, overdensity and entropy. It describes 
redshift dependence of the observed characteris- 
tics of absorbers discussed in Sec. 4, discrimi- 
nates between the random and regular variations 
of absorbers characteristics and allows one to ob- 
tain a reasonable statistical description of these 
characteristics. It also explains the observed red- 
shift variations of the linear number density of 
absorbers and links them with the fundamental 
parameters of the initial power spectrum. 

Such a complex approach allows one to test 
its self consistency through the comparison of 
several measured statistical characteristics with 
corresponding theoretical expectations. As was 
shown above, it also allows one to test the shape 
of the correlation function of the initial velocity 
field and of the initial power spectrum on unprece- 
dently small scale. 

The progress achieved demonstrates again the 
key role of the representativity and completeness 
of the observed samples for the construction of the 
physical model of absorbers and reveals a close 
connection between the observational database 
and obtained results. Some of the most interesting 
results were derived from the distribution of weak 
absorbers but unfortunately the available sample 
of such absorbers has only limited statistical sig- 
nificance. Further progress can be achieved with a 
richer sample of observed absorbers and by testing 
our approach on a more representative numerical 
simulations. 

9.1. Main results 

Main results of our analysis can be summarized 
as follows: 

1. We demonstrate that the physical model of 
absorbers proposed in Paper I and Sec. 3 is 
self consistent and agrees well with the ba- 
sic observed characteristics of the Lyman-a 



clouds. This fact indicates that the majority 
of absorbers can be considered as long-lived, 
gravitationally bound and relaxed. 

2. Using the Zel'dovich theory of gravitational 
instability we link the PDF of DM column 
density found for the observed sample of 
Lyman-a clouds with the correlation func- 
tion of the initial velocity field, what al- 
lows one to estimate the initial power spec- 
trum to unprecedently small scale. 

3. We found that at larger scales both the mea- 
sured function and the initial power spec- 
trum of perturbations do not differ by more 
than ~ 10-15% from the predictions of the 
CDM model. This result confirms the basic 
predictions of the inflationary scenario. 

4. At scales - 30 - SOO/i-^kpc we sec notice- 
able differences between the measured corre- 
lation function ^„ and the expected one for 
the CDM-like power spectrum. These differ- 
ences can be considered as an observational 
evidence in favor of a more complex infla- 
tion, which generates an excess of power at 
small scales. In turn, this excess of power ac- 
celerates the formation of first galaxies and 
causes early reheating of the Universe. 

5. Our quantitative estimates depend upon the 
limited statistics of weak absorbers and their 
poorly known survival probability. However, 

both confirmation of CDM like shape of the 
initial power spectrum or detection of its dis- 
tortions are equally important. 

These results indicate that the proposed approach 
is quite promising and the analysis should be con- 
tinued with a richer sample of observed Ly-a ab- 
sorbers. Application of this approach to simulated 
matter distribution at z = and to the sample 
of - 100 000 galaxies compiled in the SDSS DRl 
demonstrates its high efficiency for the quantita- 
tive description of strongly nonlinear matter con- 
densations. However, its application for the Ly-a 
forest at high redshifts needs to be directly tested 
with a representative numerical simulations what, 
in turn, will allow to improve the estimates of the 
initial power spectrum. 
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9.2. Physical model of absorbers 

The physical model of absorbers introduced in 
Paper I and improved in Sec. 3 links the measured 
z, b and Nhi with other physical characteristics of 
both gaseous and DM components forming the ob- 
served absorbers. In this model we consider the 
majority of absorbers as long-lived gravitation- 
ally bound and partly relaxed objects formed in 
the course of both linear and nonlinear conden- 
sation of the DM component. The existence of 
galaxies and quasars at high redshifts and the ob- 
served reheating of the Universe at redshifts z > 6 
demonstrate the importance of nonlinear processes 
at redshifts under consideration. Numerous simu- 
lations (see, e.g., Frenk 2002) also confirm the for- 
mation of high density DM filaments and sheets 
at these redshifts. 

The Zel'dovich theory of gravitational instabil- 
ity indicates the self similar character of the pro- 
cesses of DM structure formation and allows for- 
mation of numerous high density relaxed DM fila- 
ments and sheets already at high redshifts. Their 
typical sizes and their mass function depend upon 
the initial power spectrum and their mean charac- 
teristics progressively change with time. However, 
the DM structure at high redshifts is expected 
to be qualitatively similar to the one observed at 
small redshifts. The Ly-a clouds trace this DM 
structure and, as was shown in Sec. 4, the PDFs 
of their observed characteristics only weakly de- 
pend on the redshift. 

For low mass pancakes the spatial distribu- 
tion of the baryonic and DM components can be 
strongly biased due to the influence of the gaseous 
pressure (Jeans damping). The physical model of 
absorbers considered in this paper allows for this 
bias. 

Our approach is consistent with the weak red- 
shift variations of the mean observed characteris- 
tics of absorbers, namely, (6) and {Nhj), and even 
the shape of their distribution functions (Sec. 4) 
does not change with time. Moreover, the main 
statistical characteristics of absorbers discussed in 
Paper I and in Sees. 4 & 5 coincide with theoretical 
expectations. In particular, the redshift variations 
of {Q do not exceed 10% (Sec. 4), the mean value 
of C and its measured PDF coincide with the ex- 
pected ones (Sec. 5.3). The same approach allows 
one to describe the redshift variation of the linear 



number density of absorbers (Sec. 6). 

Analysis of the main characteristics of ab- 
sorbers performed in Paper I shows that the sam- 
ple of observed absorbers can be composed of 
DM pancakes with complex evolutionary histo- 
ries. They can be formed due to merging of earlier 
formed pancakes or compression of background 
matter, and both adiabatic, for (3 < 2, and shock 
compressions, for /? > 2 are equally possible. The 
great diversity of these histories is clearly seen al- 
ready from the well known broad distribution of 
absorbers in the "b — Nhi" plane (see, e.g., Hu et 
al. 1995). 

Properties of absorbers are changing even after 
formation of gravitationally bound pancakes. In 
Paper I we discussed five most important factors 
that determine absorbers evolution after forma- 
tion. They are: the transverse expansion and/or 
compression of pancakes, their merging, the radia- 
tive heating and cooling of compressed gas, and 
disruption of pancakes into a system of high den- 
sity clouds. To illustrate the influence of these 
factors and correlations between different charac- 
teristics of absorbers we roughly discriminated in 
Paper I three subpopulations of absorbers with 
various entropies and overdensities. 

Our approach allows one to discriminate be- 
tween the systematic and random variations of 
properties of absorbers. The former ones are nat- 
urally related to the progressive growth with time 
of the DM column density of absorbers, q{z), what 
can be described theoretically. The action of ran- 
dom factors cannot be satisfactory described by 
any theoretical model. However, in the frame- 
work of our approach, the joint action of all ran- 
dom factors can be described by one random func- 
tion, which can be identified with the entropy of 
the compressed gas and it is directly expressed 
through the observed parameters (Paper I). 

9.3. Reconstruction of the initial power 

spectrum 

Naturally, the formation and further evolution 
of DM pancakes is determined by the initial power 
spectriim. However, the great diversity of their 
evolutionary histories makes the reconstruction of 
this spectrum from observations of the Ly-a for- 
est difficult. In this paper we solve this problem in 
the framework of the complex statistical investiga- 
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tion of absorbers that allows one to check the self 
consistency of this approach. In spite of this, our 
results are based on some assumptions and, so, are 
model dependent. To check and to improve this 
approach we have to apply it to the high resolution 
representative simulations of formation of galaxies 
and absorbers. 

This statistical description connects the PDF 
of the DM column density of pancakes with their 
survival probability, Ws , and the correlation func- 
tion of the initial velocity field, For Gaussian 
initial perturbations and for a given Ws , this con- 
nection allows one to reconstruct the approximate 
shape of the function and to reveal its possible 
divergence from the expected one for the CDM- 
like initial power spectrum corrected for damping 
due to a finite mass of DM particles. In turn, this 
divergence can be interpreted as an indication of 
existence of special features in the initial power 
spectrum at small scales. 

In our approach we use the measured PDF of 
the DM column density of pancakes rather than 
the smoothed density field. It therefore avoids the 
complicated problem of density variations along 
the line of sight and mass conservation. It does 
not depend on the peculiar velocity of absorbers 
and is not restricted by the Nyquist limit and the 
distances between absorbers. So, this indirect ap- 
proach allows one to check the shape of the initial 
power spectrum in the real space up to unprece- 
dently large wave numbers, r] = k/ko 100. 

For larger scales, > 1 — 2hr^ Mpc, 77 < 3, our 
analysis verifies the CDM-like shape of the initial 
power spectrum and the correlation function 
How(?vc;r at small scales noticeable differences be- 
tween the CDM-like and measured spectra are de- 
tected. The reliability of these differences is prob- 
lematic (sec Sec. 7) but both confirmation and/or 
elimination of these distortions are equally impor- 
tant. 

9.4. Numerical simulations and the model 
of "fluctuating Gunn— Peterson ap- 
proximation" 

The leading role of numerical simulations in in- 
vestigation of the LSS is well known. They clarify 
many problems of structure formation and evolu- 
tion at all redshifts and, in particular, link the ob- 
served characteristics of absorbers with those that 



are not directly observable. This link allows one 
to test the physical model of absorbers and, fur- 
ther on, to use it to recover the density fields and 
the initial power spectrum from the observed Ly 
a forest. Here we shortly compare our results with 
numerical simulations. 

9.4- 1- The model of "fluctuating Gunn-Peterson 
approximation" 

One of the popular model of the observed Ly- 

a forest is the "fluctuating Gimn-Peterson ap- 
proximation" (FGPA) widely discussed during the 
last five years (see, e.g., Croft et al. 1998, 2002; 
Nusserfe Haehnelt 2000; McDonald ct al. 2000; 
Phillips et al. 2001; Zaldarriaga, Hui & Tegmark 
2001; Shaya et al. 2002). It is focused on the con- 
nection between the observed redshift variations 
of the flux emitted from quasars and the initial 
power spectrum. 

In this approach absorbers are considered as 
formed by a mildly nonlinear unshocked compres- 
sion of the background matter. It links absorp- 
tion lines with large diffuse gaseous clouds with 
moderate overdensities which are expanded by the 
residual Hubble flow. For such clouds their veloc- 
ity profile depends mainly on their size, while the 
thermal broadening is considered as minor (Croft 
et al. 2002). However, other authors consider this 
broadening as quite essential (see, e.g., Meiksin, 
Bryan & Machacek 2001; White & Croft 2003). 

The main result of the FGPA is the determi- 
nation of the power spectrum. This method does 
not use the observed hydrogen column density and 
the Doppler parameter. It begins with the analy- 
sis of the measured one dimensional spectrum of 
the flux and converts it into the power spectrum 
of the matter. The nonlinear effects are partly 
taken into account when the measured flux spec- 
trum is converted to the matter spectrum with the 
special function b{k) (Croft et al. 2002). This sim- 
ple method, in spite of its deficiencies, surprisingly 
well restores the CDM-like power spectrum down 
to scales ~ l/i~^Mpc. 

The application of the FGPA is restricted by 
the Nyquist limits. For the typical separation of 
observed absorbers (32) 

(i?,ep)«l(j^)/i-'Mpc, (49) 

the analysis allows one to reach a typical wave 
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number 

At such scales, the similarity of the observed spec- 
tra with the CDM-like one is confirmed by our 
results (44). 

Comparison of simulated and measured flux 
characteristics in McDonald et al. (2000) shows an 
excellent agreement. However, later analysis re- 
vealed systematic inconsistencies between the ob- 
served and simulated intermittent features of the 
flux that indicates the underestimation of nonlin- 
ear evolution in simulations (Pando et al. 2002). 
White & Croft (2003) note also the importance 
of the non-linearity of matter clustering for the 
reconstruction of the power spectrum. 

The FGPA has been criticized, for example, by 
Zaldarriaga, Scoccimorro & Hui (2002) who no- 
tice the limited usefulness of available simulations, 
limitations of the inversion technique and the im- 
portance of nonlinear effects. The same problems 
are also discussed in Gnedin & Hamilton (2002), 
Tegmark & Zaldarriaga (2002) and Seljak, Mc- 
Donald & Makarov (2003). 

9.4-2. Properties of simulated forest 

Usually the comparison of observed and simu- 
lated characteristics of the Ly-a forest is restricted 
by the flux emitted from quasars but always they 
are found to be quite consistent. However, the 
flux integrates the HI distribution over absorber 
and therefore it does not characterize it uniquely. 

More detailed comparison of the Doppler pa- 
rameter, 6, and the column density of neutral 
hydrogen, Nhi, performed in Meiksin, Bryan & 
Machacek (2001) reveals that no simulations pro- 
vide a statistically acceptable match to the ob- 
served b and Nhi distributions. In particular, 
they found that the median measured Doppler pa- 
rameters exceed the simulated one by as much as 
30 - 60% especially for optically thin lines and 
also the simulated and measured Nhi distribu- 
tions tend not to agree. It can be expected that 
results presented in Sec. 4.2 enhance divergences 
between the simulated and observed quantitative 
characteristics of the forest. 

As was noted in Paper I, for the sample dom- 
inated by recently formed 'young' absorbers a 



strong redshift evolution of observed characteris- 
tics and, in particular, of the Doppler parameter, 
is expected. Perhaps, this effect enhances distor- 
tions between measured and simulated Doppler 
parameters. 

The available numerical simulations are re- 
stricted by technical limitations but, in spite of 
this, they could help to solve many unclear prob- 
lems of the early structure formation. First of all, 
it is link between the gaseous and DM components 
and estimate of the degree of nonlinearity. A de- 
tailed description of the DM distribution is now 
possible with at least two effective; methods the 
Minimal Spanning Tree (Demiahski et al. 2000; 
Doroshkevich et al. 2001) and the Minkowski 
Functional (see, e.g., Schmalzing et al. 1999) 
techniques. Characteristics of DM distribution are 
extremely important for the reconstruction of the 
processes of formation of the Ly-a forest. These 
characteristics can be directly compared with ex- 
pectations of the Zel'dovich theory discussed in 
Sec. 2. 

Now it is not clear why the FGPA links the 
Ly-a forest with large diffuse gaseous clouds with 
moderate overdensities and why the high overden- 
sity DM pancakes relaxed or unrelaxed - are 
not seen in simulations. It is natural to expect 
that the formation of galaxies is accompanied by 
formation of numerous high overdensity pancakes. 
Simple estimates show that, for parameters dis- 
cussed in Sec. 4.2, the one dimensional collapse 
leading to formation of high density pancakes is 
quite probable at all redshifts under consideration. 
For {b)/bhg > Pthr shock compression dominates. 

It is also not clear what happens with the 
formed DM clouds in the course of their further 
evolution and why merging of pancakes is sup- 
pressed. It is one of the important processes re- 
sponsible for evolution of the structure which is 
clearly seen in numerous high resolution simula- 
tions performed in large boxes. This process is 
accompanied by shock compression of the gas and 
an essential growth of its entropy. In turn, vari- 
ations of the entropy explain the strong scatter 
of absorbers in the "6 — Nhi" plane (Hu et al. 
1995). This list of questions can be substantially 
extended. 

To achieve the required resolution, the simula- 
tions of the forest are usually performed in rela- 
tively small boxes with sizes Lh 

OX 

5-20/i-iMpc 
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what eliminates the large scale initial perturba- 
tions. To test how accurately such numerical sim- 
ulations reproduce observations in Fig. 10 we 
compare the velocity correlation function ^„ for 
CDM power spectrum with velocity correlation 
functions for two truncated spectra used in simula- 
tions. As is seen from this Figure, both correlation 
functions derived from truncated spectra at scales 
IvQko > 0.4 are suppressed. This in turn slows 
down the merging of clouds, decelerates the for- 
mation of richer absorbers and distorts their other 
characteristics. As was found in Sees. 4.2 and 5.3, 
for the cosmological model (1) the mean size of 
compressed absorbers is 

Hy^io.„,(_i_)\ ,50, 

what coincides with the critical value found above. 

This estimates indicate that in simulations per- 
formed with koLbox < 5 — 6 the formation of 
richer absorbers with q > {q{z)) is strongly sup- 
pressed and the weaker adiabatically compressed 
absorbers dominate. On the other hand, estimates 
discussed in Sec. 4.2 verify that at redshifts un- 
der consideration at least half of matter could be 
already accumulated by absorbers and, so, their 
evolution could be driven by the processes of merg- 
ing and shock heating. These inferences should 
be tested on simulated DM characteristics of pan- 
cakes and compared with expected ones for the 
corresponding power spectrum, for example, with 
methods described in Sec. 2 . 

9.4-3. Measurements of the power spectrum 

The FGPA picture seems to be quite different 
from our representation of absorbers as a network 
of discrete gravitationally bound pancakes with a 
variety of overdensities. Our model identifies the 
Doppler parameter with the depth of the DM po- 
tential well and the temperature of compressed gas 
rather than with the residual expansion of mat- 
ter. It assumes domination of long lived absorbers 
what agrees well with results obtained in Sec. 4.2. 
Direct estimates obtained with the pancake sizes 
and Doppler parameters specified in this section 
show that the formation of such high density ab- 
sorbers is quite probable at all redshifts under con- 
sideration. 

Both models describe quite well some of the 
observed properties of the Lyman-a forest but 



they differ in emphasizing domination of discrete 
or diffuse absorbers. In many respects they are 
complementary to each other because diffuse and 
discrete absorbers coexist and fractions of these 
types of absorbers change with redshift. Evidently 
the diffuse absorbers dominate just after reheating 
while discrete absorbers dominate at lower red- 
shifts. These differences demonstrate the complex 
character of the Lyman-a absorption clouds and 
indicate that both models should be improved af- 
ter detailed analysis of truly representative simu- 
lations of the forest. 

The proposed here method of measuring the 
power spectrum is based on the strong connection 
between the power spectrum and the mass func- 
tion of absorbers and seems to be quite promis- 
ing. It is not restricted by the Nyquist limit and 
therefore it can be extended to scales about ten 
- thirty times smaller than the method based on 
the FGPA. It is not rigidly linked to the accepted 
physical model of absorbers which can be cor- 
rected and improved in the course of further in- 
vestigations. It can be also easily generalized for 
more refined relations between the mass function 
and correlation function of velocity field. These 
corrections can be introduced after detailed in- 
vestigations of improved and truly representative 
simulations which reproduce all observed charac- 
teristics of the Ly-a forest and link them with DM 
characteristics of pancakes. 

9.5. Reheating of the Universe and the ini- 
tial power spectrum 

Recent observations of high redshift quasars 
with z > 5 (Djorgovski et al. 2001; Becker et al. 
2001; Pentericci et al. 2001; Fan et al. 2001, 2003) 
provide a clear evidence in favor of the reioniza- 
tion of the Universe at redshifts z ^ 6, when the 
volume averaged fraction of neutral hydrogen is 
found to be //f > 10~^ and the photoionization 
rate - (0.2 - 0.8) • IQ-^^s-i . These data coin- 
cide with those expected at the end of the reion- 
ization epoch which probably took place at ^; ~ 6. 
On the other hand, recent WMAP data shift the 
redshift of reionization up to z ~ 20 ± 10 (Kogut 
et al. 2003), what implies a complex ionization 
history. 

The reionization at z ~ 6 — 10 is consistent with 
the CDM-like power spectrum with the mass of 
dominant DM particles Mdm ^ 1 — 5 keV (see 
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discussion in Paper I) . Earlier reionization is com- 
patible with a larger mass of the dominant fraction 
of DM particles, Mdm > 1 MeV or with distor- 
tions of the CDM-like initial power spectrum at 
small scales, what in turn implies a complex in- 
flation and a more efficient formation of galaxies 
at high redshifts (see discussion in Spergel et al. 
2003). 

Models with complex inflation and an excess 
of power at small scales were recently discussed 
in connection with distortions of the CMB po- 
larization (Naselsky & Novikov 2002; Doroshke- 
vicli, Naselsky I., Naselsky P., Novikov 2003) and 
galaxy formation (Sommer-Larsen et al. 2003). 
This means that measurements of the initial power 
spectrum at small scale are important for many 
aspects of modern cosmology. 

The future Planck mission will provide high 
accuracy measurements of the CMB polarization 
and clarify the ionization history of the Universe. 
However, these results will give again only indi- 
rect information about the possible distortions of 
the power spectrum. This means that the mea- 
surement of such distortions with the Ly-a forest 
remains the most efficient way for investigation of 
the initial power spectrum at small scale. 

The preliminary estimates obtained above seem 
to be quite promising but must be tested with 
richer samples of observed absorbers. Such inves- 
tigations will allow one to obtain more detailed 
and reliable information and better understand- 
ing of the physical processes at high redshifts and 
the shape of the initial power spectrum. 
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Fig. 1. — The mean DM column density of ab- 
sorbers, (C), the mean Doppler-parameter, (6) (in 
km/s), the mean neutral hydrogen column den- 
sity, {logNni) and the mean absorber separations, 
{Dsep)h~^Mpc, are plotted vs. redshift z. Fits 
(31) are plotted by solid lines. 



23 




12 3 

b/<b>, Nh,/<Nh,>, D3 /<D3, > 



Fig. 2. — PDFs for the Doppler parameter, P^, 

the hydrogen cohinin density, Phi, and absorbers 
separation, Pgep, for redshift 1.7 - 2.5 (soUd hnes), 
2.5 - 3 (dot lines) and 3 - 4.5 (dashed lines). 
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Fig. 3. — Cumulative distribution functions, Wq, 
are plotted vs. q/r for Model 1 (top panel), Model 
2 (middle panel) and Model 3 (bottom panel). In 
the same panels, the best fits (38) are plotted by 
solid lines. 
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Fig. 4.— PDFs, Ng, for the Model 1 (top panel), 

the Model 2 (middle panel) and the Model 3 (bot- 
tom panel) are plotted vs. (/ {Q. Solid lines shows 
fits (14) and (15) for given by (38). 
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Fig. 5. — Correlation function of the initial veloc- 
ity field, ^v{q), for the Model 1 (top panel), Model 
2 (middle panel) and Model 3 (bottom panel). 
Dashed and solid lines show fits (38) and the func- 
tion ^CDM (10). 
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Fig. 6. — Redshift distribution of absorbers, 

(nabsiz)), for 4500 observed absorbers. Fits (16) 
for Model 1 and Model 3 are plotted by solid and 
dashed lines. 
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Fig. 7. — The PDFs of dimensionless surface den- 
sity of walls, Nq, for 3 sets of walls selected from 
two high resolution DM simulations with small 
threshold richness. Fits (14) are plotted by solid 
lines. 
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Fig. 8. Correlation functions for the initial ve- 
locity field (10) (solid line) and (38) for the Model 
1 and Model 3 (dashed and long dashed lines) to- 
gether with the same functions for spectra (47) 
and (48) (dotted lines). 
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Fig. 9. — Normalized total initial power spectrum 
for two sets of parameters (48) (solid and long 
dashed lines) and CDM-like one (47) (dashed line). 
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Fig. 10. — Normalized correlation functions of ini- 
tial velocity and density fields, &^p, vs. the di- 
mensionless point separation, qlyko, for the CDM- 
like power spectrum (thick solid line) and for trun- 
cated spectra used in simulations with 256"' cells 
and the box sizes fcoibos = 2.1 (thin solid line) 
and 3.14 (long dashed line). 
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Table 1: QSO spectra used in our analysis 
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2.6 


3.1 


256 


0302 


- 003-3'2 


3.29 


2.6 


3.1 


356 


0636 


+ 680''' 


3.17 


2.5 


3.0 


313 


1759 


+ 7544 


3.05 


2.4 


3.0 


307 


1946 


+ 766^ 


3.02 


2.4 


3.0 


461 


1347 


- 246^ 


2.63 


2.1 


2.6 


361 


1122 


-4412 


2.42 


1.9 


2.4 


353 


2217 


- 2822 


2.41 


1.9 


2.3 


262 


2233 


- 606^^ 


2.24 


1.5 


2.2 


293 


1101 


- 2542 


2.15 


1.6 


2.1 


277 


0515 


-4412 


1.72 


1.5 


1.7 


76 



1. Lu et al. (1996), 2. unpublished, courtesy of 
Dr. Kim 3. Hu et al., (1995), 4. Djorgovski et al. 
(2001) 5. Kirkman & Tytler (1997), 6. Cristiani 
& D'Odorico (2000). 



